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ABSTRACT 

Recent Suzaku X-ray observations revealed that the observed entropy profile of the intracluster medium 
(ICM) deviates significantly from the prediction of hydrodynamical simulations of galaxy clusters. In this 
work, we show that gas clumping introduces significant biases in X-ray measurements of the ICM profiles 
in the outskirts of galaxy clusters. Using hydrodynamical simulations of galaxy cluster formation in a 
concordance ACDM model, we demonstrate that gas clumping leads to an overestimate of the observed gas 
density and causes flattening of the entropy profile. Our results suggest that gas clumping must be taken into 
account when interpreting X-ray measurements of cluster outskirts. 
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1. INTRODUCTION 

In recent years, galaxy clusters have emerged as one of the 
unique and powerful laboratories for cosmology and astro- 
physics. Being the largest and most magnificent structures 
in the Universe, clusters of galaxies serve as excellent tracers 
of the growth of cosmic structures. Current generation of X- 
ray cluster surveys have provided independent confirmation 
of cosmic acceleration and significantly tighten constraints on 
the n ature of dark energy (lAllen et al.1120081 IVikhlinin et al. 



2009) and alternative theories of gravity (e.g.. ISchmidt et al 
2009). New X-ray survey missions (e.g., eROSITA) are un 



derway to improve current cosmological constraints. 

Outskirts of galaxy clusters have special importance in clus- 
ter cosmology, because they are believed to be much less sus- 
ceptible to complicated cluster astrophysics, such as radiative 
gas cooling, star formation, and energy injection from ac- 
tive galactic nuclei. Dominant physical processes in the out- 
skirts are limited to the gravity-driven collisionless dynamics 
of dark matter and hydrodynamics of the intracluster medium 
(ICM). However, until very recently, observational studies of 
the ICM have been limited to radii considerably smaller than 
the virial radius of clusters. 

Recently, X-ray observations with Suzaku have extended 
measurements of the ICM profile out t o and beyond the 
virial radius for a number of cluster s (Bautz et al.l £009} 
George et alj|200<% iReronch et al.H2009t iHoshino et aLl boiot 
Kawaharada et al. 2010; Simio nescu et al.11201 11) . While mea- 



surement uncertainties are still large, the initial results indi- 
cate that the observed ICM profiles deviate significantly from 
the predictions of hydrodynamical cluster simulations (e.g., 
George et al. 2009). In addition to testing model of structure 
formation, these new measurements are important for control- 
ling systematic uncertainties for cluster-based cosmological 
measurements. 

In the hierarchical structure formation model, clusters grow 
by accreting materials from the surrounding large-scale struc- 
ture in their outer envelope. Numerical simulations predict 
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that accretion and mergers are ubiquitous and important for 
cluster formation, and the accretion physics gives rise to in- 
ternal gas motions and inhomogeneous gas density distribu- 
tion ("dumpiness") in the ICM. Recent hydrodynamical sim- 
ulations, for example, indicate that non-thermal pressure pro- 
vided by gas motions signific antly modifies the ICM pro- 
files in the cluster outskirts (e . gjlapichino & Niemeverl2008[; 
IVazza et al.1 l2009t lLau etd] 120091) . Gas clumping, on the 
other hand, also introduces biases the X-ray measurements 
of the ICM. Mathies en et alj ( 119991) . for example, pointed out 
that gas clumping can lead to an overestimate of the ICM mass 
by ~ 10% around half of the virial radius . However, the im- 
pact of gas clumping on cluster outskirts has received little 
attention in the literature. 

In this work, we investigate effects of gas clumping on X- 
ray measurements in the outskirts of galaxy clusters. Using 
hydrodynamical cluster simulations, we show that gas clump- 
ing leads to an overestimate of gas density and causes flat- 
tening of the entropy profile in cluster outskirts (r > r2oo)- 
Our results suggest that gas clumping should be treated prop- 
erly when interpreting X-ray measurements in the envelop of 
galaxy clusters. 

The paper is organized as follows. We describe our simula- 
tions in Section|2]and define the clumping factor in Section[3] 
In Section 2] we present our results. Our main findings are 
summarized in Section 

2. SIMULATIONS 

In this work, we analyze a sample of 16 simulated groups 
and clusters of galaxies . The si mul ations presented her e are 
presented in Na gai et al.l (l2007al) and lNagaiet al. (2007b), and 
we refer the reader to these papers for more details. We briefly 
summarize relevant parameters of the simulations here. 

These simulations are performed using the Ada ptive Re- 
finement Tree (ART) jV-b ody+gas-dynamics code dKravtsovl 
119991: iKravtsov et ail 120021) . which is an Eulerian code that 
uses adaptive refine ment in space and tim e, and non-adaptive 
refinement in mass (Klypinetal. 2001) to achieve the dy- 
namic ranges to resolve the cores of halos formed in self- 
consistent cosmological simulations. The simulations as- 
sume a flat ACDM model: D, m = 1 - fi A = 0.3, D, b = 0.04286, 
h = 0.7 and a% = 0.9, where the Hubble constant is defined 
as 100/zkms" 1 Mpc" 1 , and a% is the mass variance within 
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spheres of radius 8/T 1 Mpc. The simulations were run using 
a uniform 128 3 grid with 8 levels of mesh refinement. The 
box size for CL101-CL107 is 120/T 1 Mpc comoving on a 
side and is 80/T 1 Mpc comoving for CL3-CL24. This corre- 
sponds to peak spatial resolution of ss lh~ l kpc and 5/i _1 kpc 
for the two box sizes respectively. Only the inner regions 
~ 3- 10/T 1 Mpc surrounding the cluster center were adap- 
tively refined. The dark matter (DM) particle mass in the 
region around each cluster was m p ~ 9.1 x lO 8 /z _1 M0 for 
CL101-107 and m p ~ 2.7 x IO^Mq for CL3-24, while 
other regions were simulated with a bottom mass resolution. 

To assess the effects of galaxy formation on gas clumping, 
we compare two sets of clusters simulated with the same ini- 
tial conditions but with different prescription of gas physics. 
In the first set, we performed hydrodynamical cluster sim- 
ulations without gas cooling and star formation. We refer 
this set of clusters as non-radiative (NR) clusters. In the sec- 
ond set, we turn on the physics of galaxy formation, such as 
metallicity-dependent radiative cooling, star formation, super- 
nova feedback and a uniform UV background. We refer this 
set of clusters as cooling+star form ation (CSF) clusters. We 
refer the reader to Nagai et al. (2007al) for a detailed descrip- 
tion of the gas physics implemented in these simulations. 

Our simulated sample includes 16 clusters at z = and their 
most massive progenitors at z = 0.6. The propert ies of the 
simula ted clusters at z = are given in Table 1 of lLau et all 
(2010). In this work, we use the radius r2oo in which the en- 
closed mean overdensity is 200 times the critical density at 
the redshift of the output. The cluster center is defined as the 
location of the most bounded dark matter particle. To study 
the dependence on the dynamical state of clusters, we distin- 
guish unrelaxed and relaxed clusters based on the overall mor- 
phology of the mock X-ray images discussed in Nagai et al. 
(I2007bl) . 

3. CLUMPING FACTOR 

In the ICM, X-ray photons are emitted primarily by the 
scattering of electrons off ions via thermal bremsstrahlung. 
The X-ray surface brightness is given by 

S x cx / n e n H A(T,Z)dV oc l — i / p 2 m k{T,Z)&V, (1) 

where A(T,Z) is the cooling function which depends on tem- 
perature T and metallicity Z, p e and pu are the mean molec- 
ular weight of electrons and hydrogen, m p is the proton mass, 
Pgas is the gas density, and V is the volume. Since A(T,Z) 
depends weekly on T and Z sa 0.3, the observed X-ray sur- 
face brightness profile depends primarily on the average of 
gas density squared, 

Sx(r)K(p 2 g Jr))=C(r)(p gilf ,(r)) 2 , (2) 

where we define the clumping factor as, 



Note that C = 1, if the ICM is not clumpy (i.e., a single phase 
medium characterized by one temperature and gas density 
within each radial bin). 

In X-ray analyses of galaxy clusters, it is commonly as- 
sumed that the ICM is not clumpy. By setting C = 1, the 
three dimensional gas density distribution is derived from the 
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FIG. 1 . — The radial profiles of the clumping factor for the X-ray emitting 
gas with T > 10 6 K. We show results for the sample of 16 simulated clusters 
in the CSF (top) and NR (bottom) runs. Thin lines indicate individual profiles 
for relaxed (solid) and unrelaxed (dashed) clusters. Thick lines indicate the 
median profiles of each sub-sample. 

observed X-ray surface brightness profile by inverting Eq. |2] 
However, if the ICM is clumpy, the gas density infe rred f rom 
the X-ray surface brightness is overestimated by y/C(r). In 
what follows, we investigate the clumping factor of the X-ray 
emitting ICM using hydrodynamical cluster simulations. 

4. RESULTS 
4. 1 . Clumping factor profile 

In FigureQ]we plot the gas clumping factor as a function of 
radius for the sample of 16 simulated galaxy clusters. Results 
for the CSF and NR runs are shown in the top and bottom pan- 
els respectively. Since we are primarily interested in the X-ray 
emitting ICM, we consider gas with temperature T > 10 6 K 
here. We checked that the results based on T > 10 6 K cut is 
essentially identical to the calculation that takes into account 
the energy-dependent effective area of the Suzaku XIS instru- 
ment. We thus apply the temperature cut of T > 10 6 K when 
comparing to the X-ray data, unless stated otherwise. 

Figure Q] shows very large cluster-to-cluster variations, 
ranging between C ~ 1 - 100. As we illustrate below, these 
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FIG. 2. — Median clumping factor profiles of gas with different minimum 
temperature in the CSF (left) and NR (right) runs. The dot-dashed line shows 
the median clumping factor profiles for all gas. Results for T > 10 6 K and 
T > 10 7 K are shown in dashed and solid lines, respectively. 
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FIG. 3 . — Left: the mass dependence of the median clumping factor profiles 
of gas with T > 10 6 K. Results are shown for the present-day (z = 0) clus- 
ters with M 2 oo > 10 14 /r'M (solid lines) and M 2 oo < 10 14 /r'M (dashed 
lines), respectively, in the CSF (thick) and NR (thin) runs. Right: the redshift 
dependence of the median clumping factor profiles in the CSF (thick) and NR 
(thin) runs at z = (solid lines) and z = 0.6 (dashed lines). The errorbars on 
the solid lines indicate the interquartile range. 



larger in the unrelaxed systems, due to the enhanced level of 
substructures. Since th e CSF runs suffer from the 'overcool- 
ing' problem (see e.g., iRravtsov et al.ll2005h . we argue that 
the real clumping factor should be bracketed by the results of 
our CSF and NR runs. 

To investigate the dependence on gas temperature, we com- 
pute the clumping factor profiles for several different temper- 
ature thresholds. Figure [2] shows the median clumping factor 
profiles calculated for gas with T > 10 6 K and T > 10 7 K, and 
for gas of all temperatures. For the X-ray emitting gas (e.g., 
T > 10 6 K), the clumping factor is considerably smaller in the 
CSF runs than in the NR runs. In the CSF runs, the clumping 
factor is only C < 2 at r = r2oo and C ~ 7 at r = 2r2oo, which 
is less than the clumping in the NR runs. Raising the tem- 
perature threshold further to T > 10 7 K affects the clumping 
factor of the CSF run slightly, but affects the NR runs more 
dramatically, illustrating that significant fraction of clumping 
in the NR runs is due to the gas in the temperature range of 
10 6 < T < 10 7 K. For gas of all temperatures, the clumping 
factor depends very sensitively on the cluster physics (e.g., 
cooling and star formation). In this case, the clumping factor 
is considerably larger in the CSF runs than in the NR runs. 
In the CSF runs, the clumping factor increases from C ~ 2 at 
r = 0.5r2oo to C ~ 10 at r = rwo an d C ~ 100 at r = 2r^X)- In 
the NR runs, the clumping factor is comparable to that of the 
CSF runs at r = 0.5r2oo, but they are only C ~ 2 at r = rioo and 
C ~ 10 at r = 2r2oo- 

Finally, to examine the dependence on the cluster mass and 
redshift, we compare the median clumping factor profiles for 
two mass bins and redshift bins for the gas with T > 10 6 K 
in Figure [3] In the left panel, we compare the sub-samples 
forM 20 o > 10 14 /r' Mq andM 20 o < 10 14 /i _1 M Q . In the right 
panel, we compare the results for z = and z = 0.6. In the CSF 
sample, lower mass systems show a smaller clumping fac- 
tor, because compared to massive clusters, they have a larger 
fraction of lower temperature gas that are not detectable in the 
X-ray band. In the NR sample, this mass dependence is less 
pronounced. We find little evolution in the clumping factor 
profile for all runs. 



peaks correspond to high-density, low-temperature gas asso- 
ciated with infalling groups and galaxies. The amplitude of 
these peaks are generally more pronounced in the unrelaxed 
systems. The scatter is larger in the CSF runs than in the 
NR runs, because gas condensation in infalling groups and 
satellites creates higher-density, lower-temperature gas that 
enhances the clumping factor. In the CSF runs, the clumping 
factor can exceed C > 100 in regions with prominent clumps. 

Since the scatter is highly non-Gaussian, we use the me- 
dian (instead of the mean) for a measure of a typical clump- 
ing factor in our simulated clusters. In Figure Q] we indicate 
the median clumping factor for the relaxed and unrelaxed sys- 
tems using thick-solid and thick-dashed lines respectively. In 
both CSF and NR runs, the median clumping factor is approx- 
imately unity for r < 0.5r2oo, and increases gradually toward 
r ~ f"2oo. At r = r2oo, the clumping factor is C ~ 1.3 and 2 for 
the CSF and NR runs respectively. At r > r^, the clumping 
factor increases rapidly with radius, reaching C ~ 5 and 10 
for the CSF and NR runs respectively. Note that the clumping 
factor of hot gas (T > 10 6 K) is smaller in the CSF runs than in 
the NR runs, because some of the hot gas is cooled out of the 
X-ray band (see Figure |2). The clumping factor is generally 



4.2. Effects of Gas Clumping on the ICM profiles 

In this section, we assess the effects of gas clumping on 
X-ray measurements of the ICM profiles. The top panel of 
Figure|4]shows that the observed and true median gas density 
profiles for the sample of 16 clusters at z = 0. Results are 
shown for both CSF (thick) and NR (thin) runs. The solid lines 
indicate the observed gas density profiles, derived assuming 
C(r) = 1 . The dashed lines indicate the true profiles, defined 
as the volume-averaged gas density. At all radii, the observed 
gas density profile is higher by y/C(r) compared to the true 
profile. In the CSF runs, the gas density is overestimated by 
15% at r = f2oo an d 120% at r = 2r2oo, respectively. Bias is 
larger for the NR runs, because of larger clumping factor (see 
Figures [TJ and |2}: the observed gas density is higher by 40% 
at r = r2oo ar, d 300% at r = 2r2oo, compared to the true value. 

Since the entropy of gas is defined as K = T ' jnj^' ', the 
overestimate of gas density due to clumping causes an un- 
derestimate of the observed entropy profile by C(r) 1 / 3 . This 
is shown in the bottom panel of Figure |4] in which the en- 
tropy profiles are normal ized to X200 of the self-similar model 
( Kaisen [l986l : IVbitll200 5j). We use spectroscopic-like temper- 
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FIG. 4. — Effects of gas clumping on X-ray measurements of the ICM 
profiles. Gas clumping causes the overestimate of electron number density 
profile (top) and the flattening of the entropy profile (bottom) for gas with T > 
10 6 K. Solid lines indicate observed profiles, while dotted lines indicate true 
profiles. Thick and thin lines correspond to CSF and NR runs, respectively. 



ature T s i ( Mazz otta et al.ll2004b to calculate entropyQ The true 
entropy profil e is consistent w ith the self-similar prediction 
with K oc r 1A (IVoit et al.ll2005l) . Since C(r) is radially depen- 
dent, gas clumping causes flattening of the observed entropy 
profiles at r > r^ao. 

4.3. Comparison with X-ray measurements 

In Figure|5]we compare the observed gas density, tempera- 
ture, and entropy profiles of our CS F runs to the Suzaku X-ray 
measu rements of PKS074 5-191 dGeorge etalj 120 09), Abell 
1795 dBautz et alj [2009b. Abell 1689 dKawaharada et all 
l20!ol). Abell 1413 dHoThinoetalJ l2010h. and Perseus East 
(Simionesc u et al.ll201 ll) . The temperature and entropy pro- 
files are no rmalized to T200 and ^200 from the self-similar 
model (e.g. JVoidl2005l) . 

The "observed" ICM profiles of the simulated clusters are 
generally consistent with the Suzaku measurements, with 
the exception of PKS0745-191 which deviates significantly 
from our simulated clusters as well as other Suzaku mea- 
surem ents. Since the sp ectroscopic-like temperature pro- 
files dMazzotta et ail 120041) are in good agreement (see also 
iBurns et al.l l2010j) ! the discrepancy in the entropy profile 
arises primarily from the difference in the gas density profile. 

3 We only consider X-ray emitting gas with T > 10 6 K when computing 
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FIG. 5. — Comparisons of simulated clusters and Suzaku X-ray measure- 
ments. From top to bottom, we show the radial profiles for the electron num- 
ber density, spectroscopic-like temperature, and entropy. For the simulated 
clusters, we show the median profiles of the CSF runs. Dashed lines represent 
the true profiles. Solid lines represent the observed profiles, with thin-solid 
lines indicating th e interquartile range. Black points are S uzaku observations 
of PK S0745-191 (George et al. 2009) (ci rcles) A1689 )Kawaharada"etaI] 
120101) (triangles) A14 13 <Hoshino etal] 12010) (stars), and Perseus East 
I Simionescu et al. 201 1) (squares). The black dot-dashed lines are the best-fit 
profile of A1795 I Bautz et al. 2009) and the dotted black lines indicate ±20% 
uncertainty around the best-fit model. 



While the flattening of the median profile of our simu- 
lated clusters is not as significant as the observed profile of 
PKS0745-191, we note that the simulated clusters show sig- 
nificant cluster-to-cluster variations, owing to the large varia- 
tions in the clumping factor shown in Figure [U. In fact, sev- 
eral simulated clusters exhibit the level of flattening similar 
to that of PKS0745-191 in their outskirts. For example, as 
shown in Figure Q] one of the relaxed clusters, classified as 
relaxed based on its morphology in the inner regions, exhibits 
a large clumping factor of C ~ 100 at r = ^200- 
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However, we caution that these measurements are still quite 
uncertain. For example, it is puzzling that the gas density 
of PKS0745-191 is considerably larger than the other Suzaku 
measurements at all radii. Observations of more clusters and 
better understanding of systematic uncertainties are required 
before making robust conclusions. Independent measure- 
ments with Chandra, ROSAT, Swift-XRT, and XMM-Newton, 
can shed light o n systematic uncertainties in the current mea - 
surements (e.g., Ett ori & Molendill2010HMoretfi et ail 2010). 



5. CONCLUSIONS AND DISCUSSIONS 

In this work, we show that gas clumping can cause signif- 
icant biases in X-ray measurements of ICM profiles. Using 
hydrodynamical cluster simulations, we show that gas clump- 
ing leads to an overestimate of the gas density and causes flat- 
tening of the entropy profile in the outskirts of galaxy clusters 
(r > r 2 oo)- 

By comparing results of simulations with different input 
cluster physics, we find that the magnitude of gas clump- 
ing depends sensitively on the input cluster physics (e.g., 
gas cooling and star formation). For the X-ray emitting gas 
(T > 10 6 K), the amount of gas clumping is smaller in clus- 
ters with radiative cooling and star formation (CSF) than those 
with simple non-radiative (NR) gas physics, because cooling 
removes some of the high-density, X-ray emitting gas out of 
the detectable temperature range. Our hydrodynamical simu- 
lations show that C ~ 1 .3 and 2 at r = ^200 and C ~ 5 and 10 at 
r = 2r2oo for the CSF and NR clusters, respectively. Since our 
CSF runs suffer from the overcooling problem, we argue that 
the real clumping factor should be bracketed by the results of 
our CSF and NR runs. We also checked that our results are es- 
sentially unchanged for the simulations for the WMAP-5 cos- 
mology: n M = 1 - fl A = 0.27, fl B = 0.0469, h = 0.7, cr 8 = 0.82. 

The simulations also show that there is a large variation of 
the gas clumping factor among the cluster populations. The 
gas clumping is more pronounced in dynamically active sys- 
tems, while it depends mildly on the cluster mass and redshift. 



Given the scatter, our results are consistent with the current 
observational constraints (e.g., Sim ionescu et al. 2011). 

While current Suzaku measurements are still uncertain, the 
present work suggests that gas clumping is important for re- 
ducing the tension between the recent Suzaku observations 
and the theoretical predictions of the ACDM model. We 
therefore caution that the effects of gas clumping must be 
taken into account when interpreting X-ray measurements of 
cluster outskirts. Results of this work will also be relevant for 
interpreting future X-ray missions, e.g., eROSITA and WFXT, 
which will produce large statistical samples of ICM measure- 
ments out to large cluster-centric radii. 

Further studies of gas clumping will be important for testing 
the predictions of hydrodynamical simulations and for using 
galaxy clusters as precision cosmological probes. Analyses 
of the nature of gas clumps (e.g., size distribution, filling fac- 
tor, and angular dependence) using hydrodynamical simula- 
tions will be useful for designing observational strategies to 
constrain gas clumping. For example, gas clumping may be 
enhanced along the directions of filaments where the cluster 
accretes clumpy and diffuse materials, and it might help ex- 
plain the anisotr opic distribution o f the ICM seen in recent 
simulations (e.g. .[Vazza et aLl uOlO) and in observations (e.g., 
iKawaharada etaiT T20 10). Since the Sunyaev-Zel'dovich ef- 
fect (SZE) signal is less sensitive to gas clumping, compar- 
isons of X-ray and SZE measurements will shed light on the 
physics of cluster outskirts. 
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